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Abstract. A symmetry-breaking phase transition in the early universe could have led to 
the formation of cosmic defects. Because these defects dynamically excite not only scalar 
and tensor type cosmological perturbations but also vector type ones, they may serve as a 
source of primordial magnetic fields. In this study, we calculate the time evolution and the 
spectrum of magnetic fields that are generated by a type of cosmic defects, called global 
textures, using the non-linear sigma (NLSM) model. Based on the standard cosmological 
perturbation theory, we show, both analytically and numerically, that a vector-mode relative 
velocity between photon and baryon fluids is induced by textures, which inevitably leads to 
the generation of magnetic fields over a wide range of scales. We find that the amplitude of 
the magnetic fields is given by H ~ 10“®((1 -|- (u/mpi)^ (fc/Mpc“^)^'^ /VN Gauss 

in the radiation dominated era for A: < 1 Mpc“^, with v being the vacuum expectation value 
of the 0(N) symmetric scalar fields. By extrapolating our numerical result toward smaller 
scales, we expect that B ~ 10 “^^ ((1 -|- 2 :)/ 1000 )~^^^ {v/rripif' (/c/Mpc ”^)/\fN Gauss on 
scales of A: > 1 Mpc“^ at redshift z > 1100 . This might be a seed of the magnetic fields 
observed on large scales today. 
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1 Introduction 

The standard big-bang cosmological model has been established by a variety of astronomical 
observations, such as the dimming of distant supernovae, the cosmic microwave background 
(CMB), the large-scale structure of the universe, and so on. In this big-bang paradigm, it 
is naturally expected that the universe experienced a number of phase transitions at the 
early stages of cosmic history, because of cooling from a very hot initial state caused by its 
adiabatic expansion. As a result of these phase transitions, topological defects may have been 
formed and may have acted as a source of large scale structure in the universe. Among others, 
topological defects resulting from a phase transition that breaks a global 0{N) symmetry are 
called 0{N) global defects. For example, 0{N) symmetric scalar fields can yield domain 
walls (A = 1), monopoles {N = 3), and textures (A > 4) when the scalar fields break the 
global symmetry. 

There exist many studies regarding the phenomenological aspects of topological de¬ 
fects, which include the generation of CMB temperature and polarization anisotropies [1-4], 
gravitational waves (GWs)[5-9], and cosmic rays [10, 11], among others (for a review of struc¬ 
ture formation with topological defects, see [12]). As a rule of thumb, the amplitude of 
the fluctuations induced by topological defects, such as CMB anisotropies, is of the order 
AT/T ~ AttGv^ , where v is the vacuum expectation value (VEV) of the scalar fields, G is the 
Newton constant, and T and AT are the CMB temperature and its fluctuation, respectively. 
Therefore, the recent CMB measurement by the Planck satellite has placed limits on the 
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energy scale of the topological defects [13]. The actual limits depend on the detailed models, 
for instance, < 4.2 x 10“^ [14] for cosmic strings, < 3.2 X 10“^ for Abelian-Higgs 
cosmic strings, Gv"^ < 1.5 X 10“^ for Nambu-Goto strings, and Gu^ < 1.1 X 10“® for semi-local 
strings and global textures [13]. 

In this paper, we pay particular attention to global textures with S> 4, which can 
be well-approximated by self-ordering scalar fields that follow the non-linear sigma model 
(NLSM). The NLSM can describe the evolution of global 0{N) symmetric scalar fields with an 
accuracy up to corrections on the order of 1 /N. The NLSM has attracted much attention since 
the discovery of CMB B-mode polarizations by the BICEP2 experiment [15] because textures 
following the NLSM can be a source of the scale-invariant GWs [16] [7] [5], just as inflation 
in the early universe can produce the scale-invariant GWs. To observationally distinguish 
between GWs originating from inflation and from textures, one should consider observables 
that reflect the time evolution of the GWs. The GWs from inflation are frozen on super¬ 
horizon scales at first, and decay with oscillations after the horizon crossing. The GWs from 
textures, on the other hand, are generated inside the horizon, and decay with oscillations 
after the scalar fields that source the GWs decay away as the universe expands. In [17], 
the authors calculate the CMB temperature and polarization anisotropies in the NLSM and 
find that the shapes of the correlation functions of the CMB anisotropies in the NLSM are 
different from the corresponding ones from inflation. Therefore, detailed observations of 
CMB anisotropies can distinguish between GWs from the two different origins. In fact, the 
recent B-mode measurement by BICEP2 places an upper bound on the VEV in the NLSM at 
u < 9 X 10“^G“^/^, and the GWs from the NLSM are shown to be slightly disfavored by the 
data compared with those from inflation [18] [19], while contamination by dust in the BICEP2 
data has to be re-analyzed with PLANCK data. 

In this paper, we investigate yet another route to probe the NLSM: the generation 
of magnetic fields. Because of the non-linear nature of the NLSM, scalar fields following 
the NLSM inevitably induce vector-mode perturbations as well as scalar and tensor-mode 
ones. The relative vector-mode velocity between photon and baryon fluids induces rotation 
in electric fields, leading to the generation of magnetic fields [20]. Recent discoveries of large- 
scale magnetic fields in void regions [21] [22] as well as magnetic fields at high redshifts [23] 
make the investigation more interesting because such magnetic fields may be of primordial 
origin in the early universe. Therefore, one of the aims of this paper is to derive the spectrum 
of magnetic fields generated by vector perturbations in the NLSM within the observational 
limits of CMB anisotropies. 

This paper is organized as follows. In the next section, we review the NLSM, in which 
Wcomponent scalar fields act as a source of cosmological perturbations. In section 3, we 
derive the power spectrum of magnetic fields in the NLSM, both numerically and analytically, 
using the tight coupling approximation. We will see that in order to obtain a reliable result, 
we should expand the equations up to the third order in the tight coupling approximation 
between photon and baryon fluids. We discuss the result and give an analytic interpretation 
of the spectrum of the magnetic fields in section 4, followed by our conclusion in section 5. 
Throughout this paper, we fix the cosmological parameters to h = 0.7, = 0.0226, = 

0.112, and = 3.046, where Hq = lOO/i km/s/Mpc is the Hubble constant, H;, and He are 
the density parameters of baryonic and cold dark matter, respectively, and Ny is the effective 
number of massless neutrinos. Those parameter values are consistent with the Planck results, 
and correspond to the AGDM model [24]. 
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2 Non-linear sigma model (NLSM) 


We consider a model of N scalar fields with a global 0{N) symmetry, which undergoes 
symmetry breaking in the early universe. The Lagrangian of the model is given by 

£ =( 2 . 1 ) 

Here = {4>i,(l)2, ■■■, 4 >n) is an array of N real scalar fields, v is the vacuum expectation value 
(VEV) of the scalar fields after symmetry breaking, A is the self-coupling constant, and Lt ~ 
is the thermal term of the Lagrangian. Deep in the radiation-dominated universe, 
the thermal term is dominant in the Lagrangian of the scalar helds. As the temperature of 
the universe goes down, this thermal term gets smaller. When the term becomes negligible 
a spontaneous symmetry breaking occurs. At energy scales well below v, the field values 
are confined on the N — 1 dimensional sphere in the N dimensional field space so that 
? 7 )<h(x, r/) = v'^. Under this condition, in the large-limit, the equation of motion for 
the scalar fields can be derived as 

V^V^/3a + (V^/3b) • {W/3b) 13a = 0, (2.2) 

where the indices a and b run over 1,...,A^ and summations with respect to deprecated 
indices are implicit. This is the equation of motion for the scalar fields in the NLSM. If 
N <3, topological defects such as domain walls, cosmic strings and monopoles are produced, 
where the scalar fields can restore the 0{N) symmetry and possess the higher energy density. 

In this paper we consider cases with > 4, where non-topological defects, or textures, 
form. Let us consider the flat Friedmann-Lemaitre-Robertson-Walker universe with the metric 

ds^ = = a^{ini){—drf‘ + ddP') , (2.3) 

where 0 ( 77 ) is the cosmic scale factor and r/ is the conformal time. In the large-A^ limit, making 
an ansatz that ((V^/3a) • (V^/3a)) = we can obtain the analytic solution as [7] [5] 

I3a{k, rf) = ^ V*)l3a{k, 77 *), (2.4) 

where u = dlna/dln 77 -|- 1, = 4r(2;z — l/2)T{v — l/2)/3r(z^ — 1) and Tq = 3v — 9/4. In 

eq. (2.4), 77 * is the time of phase transition and /3a(^)^*) is the initial value of scalar fields, 
whose two-point correlation function can be given as 

(^l3aiki,r]^)Pl{k2,V*)') = ^ab {2TTf6{ki - ^ 2 )- (2.5) 

The above relation is only valid for k <C I/ 77 *, which follows, in the large-A^ limit, from the 
fact that the scalar fields take their VEV independently in each horizon at 77 *.^ We also note 
that I3a{k,r]^) is Gaussian at these scales. 

We denote a correlation function of f3a as a product of the transfer function and the 
initial amplitude as 

(^f3a{h,v)Pbik2,v)) = F'{ki,k2,rj,r]^) (^(3aiki,v*)l3h{h,V*)) , ( 2 - 6 ) 

^ On small scales k > l/r],, initial scalar fields become homogenous and correlation function of eq. (2.5) 
should vanish. We however note that our argument does not depend on the correlation function on these 
scales. 
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where F{ki,k 2 ,'r],'t]*) = f{ki,r],r]^)f[k 2 ,r],r]^). From eq.(2.4) and eq.(2.5), we can see that 
the dependence of eq.(2.6) on rj^, is canceled out. Therefore we omit r/* from the equations 
hereafter and write as F{ki,k 2 ,Tj,r]^,) = F{ki,k 2 ,Tj) and = Pa{k). Finally the 

energy-momentum tensor of the scalar fields is given by 


rp4, 

^ flU 


= V 


di^PaduPa 


^giMudxPad Pa 


(2.7) 


3 Magnetic fields 

In this section, we investigate generation of seeds of large scale magnetic fields from the self¬ 
ordering scalar fields which follow the NLSM. These scalar fields can induce cosmological 
vector-mode perturbations and eventually produce magnetic fields. 


3.1 Vector mode perturbations and their evolution equations 

We begin by reviewing the basic linear perturbation theory and define the vector mode. Let 
us consider the perturbed metric around the flat FRW one in the synchronous gauge as 

ds"^ = g^iudx^dx^ = a^{r]){—dif‘ + {5ij -|- hij)dx^ dx^), (3.1) 


where hij is the metric perturbation. In Fourier space, the vector part of hij can be expressed 
as 


hij 


ikihj + ikjhf 

71 


(3.2) 


Here hf is a divergenceless vector and can be rewritten using the vector basis e\^\k) as 


hY = Xh\Y^e\^\k). 

A=± 


(3.3) 


Combining eq. (3.2) and eq. (3.3), we can denote hij directly as 

= Y ^ Y (3.4) 

A=± 

where is the vector projection tensor. Using this projection tensor, we can derive the 
vector-mode perturbation equation for a = hy/k as 

+ 2'Ha^^'^ = SvrGa^n^^) /k, (3.5) 

where H^^) = is the total anisotropic stress in the vector mode,"!^ = a/a is conformal 

bubble, and a dot denotes a conformal time derivative. Hereafter we shall omit the superscript 
(A) for the purpose of presentation. The total energy-momentum tensor consists of two parts: 
one is from the ordinary matter and radiation and the other is from the scalar fields given by 
eq. (2.7). Anisotropic stress of the scalar fields H'^ = Tf-Oij can be calculated as 

a^n‘^(fc,r?) = y J -^^-^J^VT^^[k-2q^i\qF{q,p,r])pa{p)Pa{q){27rf6{k-p-q), 

(3.6) 
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where ^ = k-q and p = \Jk‘^ — 2kqp + For the expression of anisotropic stress of ordinary 
matter and radiation, we refer to, e.g., ref. [25]. Let us define the transfer function for the 
anisotropic stress of the scalar fields as 


z^7r‘^(A:, q, p, rf) = yx/l - [k - 2qp\ q F{q,p, ??), 


(3.7) 


and the transfer function for a as 

dvr 


dvr n _ 

a{k,q,p,r]) = / dr]' a^{r]')y'l - [k - 2qfj]q F{q,p,r]')/k. 

a^{r]) J 


(3.8) 


The fluid equation for baryon in the vector mode is given by 


Vb + Tivk 


^P'y I \ 

-^aneCFTiv^ -Vb), 

^Pb 


(3.9) 


where and pb are the energy densities of photon and baryon fluids, respectively, R = 
Apr^/Zpb, ctt is the Thomson scattering cross section, rig is the electron number density, and 
arieCTT = f is the opacity of the Thomson scattering. On the other hand, the vector-mode 
Boltzmann equation of photon fluid expanded in terms of multipole momenta is given by 


-h -A;n..y = -f{v^ - Vb), 


n.y -I- -kis - -kv^ = 
0 0 


ii + k 


i 


1 + 2 

21 +1 \TTi 


Ii+i — h-i I — 


9 ^ 9 ^ \ 8, 

' To"-'- + 5*^"' 

—fli (for I >3), 


(3.10) 

(3.11) 

(3.12) 


for intensity and 

• (/-b 3)(^-b 2)/(/- 1) ^ I 

(/ +1)3(2/+ 1) ^ 

Bi + - 


(/ +1)3(2/ + !) 


2 / + 1 
/ 

2 / + 1 


kEi_i = 
kBi_i = 


-f ( £, - d^*,) + (3.13) 

kEi, (3. Id) 


/(/ + 1 ) 


for polarization. Here and n..y are the velocity and anisotropic stress of photons, respec¬ 
tively, Ii is the /-th order moment of photons’ distribution, and Ei and Bi are the photons’ 
polarization moments and ^ = 3 / 2 /d + 9 E 2 I 2 [26]. 

As we shall show in section 3.2, the relative velocity between the photon and baryon 
fluids plays the key role in generation of magnetic fields. Since the strength of the coupling 
between photon and baryon velocities significantly changes before and after recombination, 
evolution of the relative velocity and hence the magnetic fields qualitatively differs between 
these two epochs. Before recombination, the tight-coupling approximation allows us to solve 
the system of equations partially, which we shall see shortly, and helps us to interpret nu¬ 
merical results. On the other hand, after recombination, the system of equations is solved 
completely numerically. 
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3.1.1 Tight-coupling approximation 

In the early universe, photon and baryon fluids are tightly coupled because the opacity of 
the Thomson scattering f is very large. Therefore we can expand the perturbation equations 
in section 3.1 in terms of the tight-coupling parameter k/f <C 1. This is called the tight 
coupling approximation (TCA). In ordinary analyses without external sources such as the 
NLSM, the first order solution for anisotropic stress of photons Ilg^^ and the second order 
solution for the relative velocity between the photon and baryon fluids 
were used [20]. However, when there exist NLSM scalar fields as an external source in the 
linearized Boltzmann system, we find that one should consider the TCA up to the third order 
terms proportional to cr, as discussed below. 

In the tight coupling expansion, the baryon velocity relative to the photon velocity is 
expanded using the tight coupling parameter, i.e., v^—Vb = + ..., where 

and are proportional to {k/f) and (k/f)'^, respectively. The tight coupling solutions of 
n..y and 6v up to the second order are given by 


n(b 

7 

n(2) 

7 


15 

15 






(3.15) 




n 


,( 0 ) 


fj {1 + R)k ^ 


n 


f J {l + R)k ^ 
2 


,(i) 


15 


1 + i?^^ ^ ^ 




UR n 


f J (l + ii)2A;2 \ l + R n 


+ T7+^ + 


7,(0) 

,(o) 


(3.16) 


(3.17) 


It is imporant to note that is not necessarily smaller than 5n(^) in the NLSM. This 

is because, in the NLSM, the metric perturbation a is always much larger than the fluid 
perturbation variables such as (see figure 1 and discussion in section 4), and it sometimes 
happens that the first order solution proportional to is smaller than the second order 
solution proportional to a i.e. {k/f)v^ < {k/f)‘^a. We can see directly this relation from 
(3.10), which implies 


■ / 


kr] [ ^ ] a. 


(3.18) 


Therefore, the condition that {k/f)v^ < {k/f)‘^a is satisfied at least on super-horizon scales, 
and the slip term is dominated by the second order terms in the tight coupling approximation. 

In fact, in the numerical calculations, we must eventually switch to evaluate the slip 
term directly from and Vb because the tight coupling approximation breaks down in late 
times. Because the slip term is dominated by the second order term (5n(^) in early times, we 

must keep solving Vj and Vb accurate enough up to the second order in the tight coupling 
(2) (2) 

approximation, i.e., v/ ' and \ In terms of the tight coupling approximation, the evolution 
equation of v^\ for instance, is given by 


4^) = —— T(fn(^) 


(3.19) 
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This is why we need to consider TCA up to the third order. The slip term at the third order 
is given by 




n 


fj {1 + R)k 


- 


15 


1 


1 + i? 




Tin 


( 1 ) 


(l + i?)2A:2 




nn 


+ TT V 


15 

2 

4^ 


H 


fJ {l + R^k 

3 


a 


(0) 


1 


fJ (l + i?)2 


(23 + llii)^a(°) - (17 + lli?)d(°) 


(3.20) 


l + R 


.(3.21) 


Here we show only the terms proportional to a. The condition <C is always valid 
in the tight coupling regime. 



Figure 1. Time-evolutions of the transfer functions. Here a{k, q, ii, rf) (red line), v^{k, q, /i, r]) (green), 
Vb{k^ q, /r, rj) (blue) and Sv{k^ q, /r, rj) (magenta) are plotted as functions of the scale factor. We assume 
k = q = 10“^Mpc“^ and, /i = 0. We can see that the condition cr ^ is almost always satisfied. 

3.2 Magnetic field generation 

We consider generation of magnetic fields originated from the relative velocity between the 
photon and baryon fluids, 5v = — Vb- Well before recombination, due to the frequent 

Thomson scattering of photons off electrons, electrons are separated with protons, and move 
together with photons. For protons to catch up with electrons, electric fields are induced and 
rotation of the induced electric fields generates magnetic fields via Maxwell equations. 

The equation for the generation of magnetic fields is given by [27] 

- Vbj), (3.22) 

a d ?7 3e 
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where e is the elementary charge and is the Levi-Civita tensor. The appearance of the 
rotation of 6 v in eq. (3.22) clearly shows that only the vector-mode part of Sv can contribute 
to magnetic fields. By integrating eq. (3.22) in Fourier space, we obtain 




Jo 


(3.23) 


Next we take an ensemble average of this expression over the initial configuration of the 
NLSM scalar fields j3a{k)- The ensemble average of the relative velocity can be calculated 
using the transfer function 6 v{k,q, fi,ri) and the NLSM’s initial power spectrum defined 
in appendix A.l as 




N 


Svj{k,r]')6vt{k',T]'')j = Pji{k)^{27rfd{k - k') 

X dqq^ dfi 6v{k,q, p,ri')Sv{k,q, (3.24) 


Pji{k) = 5ji - kjki- 

The correlation function of magnetic fields is then obtained as 

(^B\k,r,)Bt{k',p)) = i2nfSB{k,p)6{k-k'), 

where 


(3.25) 


(3.26) 


drj' a'^{r]')p^{r]')6v{k,q,p,r]') 


= jdq^ j d,, 

(3.27) 

We calculate eq. (3.27) numerically and the power spectra of magnetic fields at several red- 
shifts are depicted in figure 2. 


4 Analytical Understanding 

Let us try to understand the results obtained in the previous section analytically. In our 
numerical calculations, we consider scalar fields following the NLSM as the only source of 
vector mode cosmological perturbations. Therefore we assume that there are no vector mode 
perturbations at r/ —?■ 0. In this setup, vorticity a evolves first with the scalar fields as an 
external source (eq. (3.6)), and it induces photons’ anisotropic stress 11.^. Then 11.^ leads 
to the photon velocity and it propagates to the baryon velocity Vi,. Because the induced 
anisotropic stress 11.-^ is suppressed by a factor of the tight coupling parameter k/f due to 
the frequent Thomson scattering (see eq. (3.15)), the velocities induced from the anisotropic 
stress are much smaller than a, in other words, the condition that a ^ ^ Vb is valid at 

least in the tight coupling era. Therefore, we can find that the dominant term in eq. (3.17) is 


6 v{k,q,p,p) ^ 

lb \t 


1 


l + R 


a{k,q,p,r]). 


(4.1) 








Figure 2. Power spectra of magnetic fields at z = 2000 (red solid line), z = 5000 (blue solid line) and 
z = 10000 (green solid line) from scalar fields in the NLSM. Here we take the NLSM parameters as 
/N = 10“^^mpj. The black dotted line is the approximate amplitude which is given by eq. (4.18). 


To calculate a{k,q, we introduce approximations for the Bessel function in eq. (2.4) 

which are given by 


JyiX^ 


2'^T(v+l) 

^"v/Scos (x- ^vr) 


(for X ^ 1) 
(for X 1:^ 1) 


(4.2) 


Using these approximations, we can calculate super-horizon {krj <C 1) and sub-horizon [kr] S> 
1 ) solutions in the radiation and matter dominated eras as we shall show below. 

4.1 Super-horizon 


On super-horizon scales, the wavenumber of fluctuations k is smaller than the inverse of the 
horizon scale, i.e., kr] <C 1. However the wavenumber q in eq. (3.27) which comes from the 
convolution integral does not necessarily satisfy qr] <C 1. We know from eq. (3.7) and eq. 
(4.2) that a{k,q, fj,,!]) as a function of q decays as <7 oc ^ ^ rj~^, and we could 

expect that the bulk of the (/-integration comes from the range of (/ < rj~^. In this case we 
can express eq. (3.27) as 

a^-^Ssik) oc [iqa~‘^5v{k, q = 1/r/, fi, t])] ^ . (4.3) 

For k q ^ p ^ 1/t/, we can find the A:-dependence of a from eq. (3.8) as, 

cToc^. (4.4) 


Then, substituting eq. (4.1) and eq. (4.4) into eq. (4.3), we obtain 


r.3 

— T^Ssik) oc k^ 
27r^ 


e 

k 


oc k‘. 


(4.5) 
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We can see this power law tail on large scales in figure 2. The spectral index is same as that 
of the magnetic fields generated from second order density perturbations [27], but slightly 

different from the one obtained in the Einstein-aether gravity model, where ^oc k‘^ 
[28]. 


4.2 Sub-horizon 

On sub-horizon scales, i.e. fer/ ;§> 1, the situation changes in a more complicated way. Because 
of the condition that prj = — 2 kqp + q^rj 3> 1, we must take special care of rapid 

oscillations of the Bessel functions in eqs.(3.8) and (3.27). In order to manipulate the equations 
analytically, we divide the interval of integration with wavenumber q into three regions: (i) 
k ^ q ^ rj~^\ q is on sub-horizon scale, (ii) q <C q is on super-horizon scale, and (iii) 

r]~^ < q < q is nearly on the horizon scale (with a being 0(1) constant). Considering 

contributions from each interval, we estimate the power spectrum of magnetic fields. 

4.2.1 case (i) k ^ q':$> r]~^ 

In this case, because the conditions that prj ^ 1 and qrj ^ 1 are satisfied from the conser¬ 
vation of the momentum, the source function of the vorticity a has decayed away. Without 
sources, the vorticity also decays and therefore the contribution from this part is negligible. 


4.2.2 case (ii) q r] ^ 

In this case, p ^ k and the source of the vorticity is growing. Using the approximation of 
eq. (4.2), and assuming the radiation dominated era {y = 2), we can evaluate eq. (3.8) as 




O ^ j\n" 


oc 7]' ^k ^q{kp'YJ^{kr]'). 


(4.6) 


Here we use the fact that 1 — 2qp/k ~ 1 and the formula f dx x^~^^Jn{x) = x^~^^ Jn+i{x), 
and ignore the factor — p?. Using the above formula again and eq. (4.1), we obtain 

dp' p^a^Sv oc k~^q{kp)'^Ji{kp). (4.7) 

Substituting eq. (4.7) into eq. (3.27) and using eq. (4.2), the spectrum can be written as 



rl/r) 

-^Ssik) J dq q‘^ [k~'^{kp)'^J4{kp)Y 


oc k^ 


(4.8) 


4.2.3 case (iii) p ^ < q < ap ^ 

In this range, we need to treat the source carefully. First, assuming the radiation dominated 
era, let us divide the p' integral in eq. (3.27) as 


rv j-ri rq ^ p 

/ a^{p')p^{p')Sv dp'(X p'a{p')dp'= p'a^,^g-i{p')dp'+ p'a^,^g-i{p')dp' . 

Jo Jo Jo Jq-^ 

(4.9) 
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The integrand of the first term in the above equation, is given by eq. (4.6), and 

the integration leads to the term proportional to k~^q{k/q)'^Ji{k/q). That of the second term 
is calculated as 




+ 

oc 


dr] rj — 2 q]i/k) 


f 

■ f dr]"r]''^q{l — 2q]i/k) 
Ja-^ 


Jv{kr]" 


{kr]"Y 

Ju{kr]'') Ju{qd") 

{kri"Y {dY'Y 


Y ^[q{k/qfJYk/q) + k^q ^{kY){h{kY)J2{qY)+C>{q/k))] (4.10) 


where we used Ju{qY')/{qY'Y — QV" omitted the constant factor of 0(1)- 

Then the integration of 5v can be calculated, by integrating by parts, as 

[ dq' p^a^dv oc [k~^q{k/qYjA{k/q) + k~^q{k/qY JYk/q){q^ - q~Y 

Jo 

+ k~^q~^{kqYj4{kq)J2{qq)] • (4.11) 

Substituting the above equation into eq.(3.27) and ignoring the cross terms, we obtain, 
k3 _ fci/v 


^Ssik) oc dq q^ [k q^{k/qf Jl{k/q) 


+ k ^q^{k/qfjl{k/q){Y^-q ‘^f + k ^q '^{kqYJl{kq)j‘^{qq)] . (4.12) 


Integrating with g, we find that the first, second, and the third terms in the above equation 
give terms oc k^, oc k^, and oc k^, respectively. Taking these terms together, we can find the 
k dependence of the magnetic field spectrum as 

Y-^Ssik) CK k[l + 0{l/kq)]. (4:T3) 


Therefore, in the radiation dominated era {v = 2), we find a'^k^SB{k)/2TT‘^ oc k, which is 
confirmed in our numerical calculation. Finally, by reading off the numerical amplitude from 
the result of our numerical calculation we find the power spectrum of magnetic field on sub¬ 
horizon scales as 


10-^^ 1 / V \^ ( k 

(1-I-z)3 iv \^10“3?7ipjy VMpc“^ 


1 


[G^] 


(4.14) 


Similarly, we find Yk^SB{k)/2TT‘^ oc k~^ in the matter dominated era {v = 3). 

In figure 2, we find that the spectrum shows k^SB oc k on small scales (say. A: > 0.1 
Mpc“^) in the radiation dominated era. Ai z = 2000 (the red solid line), the spectrum on 
sub-horizon scales shows the k dependence between the fully radiation dominated (oc k) and 
matter dominated ones (oc k~Y- On much smaller scales {k ^ 1 Mpc“^) and in the matter 
dominated era, we expect that the spectrum of magnetic fields should be proportional to k 
because on those scales the source of vector perturbations has already decayed away and the 
magnetic fields just decay adiabatically after their creation deep in the radiation dominated 
era. From the fact that the generation mechanism is based on the mass difference between 
positively charged particles (protons) and negatively charged particles (electrons) and the 
small velocity slip between these particles, we expect that the spectrum continues up to the 
horizon scale at the epoch of annihilation k 10^*^ Mpc ^ and a cutoff at that scale [20]. 
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4.3 Approximation at super-horizon scale 

On super-horizon scales and in the radiation dominated era, we can estimate not only the 
shape of the spectrum but also the amplitude of magnetic fields approximately. On super¬ 
horizon scales the power spectrum of magnetic fields is given by 


^3 ls / Arrrr\ ^ \ n 

^ ^ —j -^k'^ dq J dr]'a^{r]')p^{r]')6v{k,q, iA,r]') 


k^ 

vr 


r-l/f? 


27r2 \ 3e 

Substituting eq. (4.2) and eq. (3.7) to eq. (3.5) we get, 

ttA 


a{k,q,p,r]) ~ -/r^ ( — ) 


48 


\mpij 


k 


(4.15) 


(4.16) 


where we assume k q ^ l/r]. Using the above expression, we can write the velocity slip as 


6v{k,q,p,r]) 


ttA, 


180 


r,Vk. 


Substituting eq.(4.17) into eq.(4.15), we obtain 


a^{z)^SB{k,z) 


10 


-3 


1 


(alA 


\mpij 


k 


[G^ 


(4.17) 


(4.18) 


27r2 ’ {1-\-zY N \mpi J \Mpc“^^ 

in unit of Gauss at redshift 2 ; ^ Zgq ~ 3300 [24]. This analytic power is plotted in figure2 to 
make a comparison with numerical results. 


4.4 After recombination 


In the above sections, we discussed magnetic field generation in the era when the tight coupling 
approximation is valid. Here we consider the era after the tight coupling approximation 
breaks down. In particular, we estimate when the magnetic fields become source-free on 
super-horizon scales, by considering the time evolution of the source function of magnetic 
fields S{r]) = pyu'^dvq^^'^, which satisfies B 2 oc [JdpS]^. 

On super-horizon scales {k <C T-L), the scalar fields with wavenumbers g ~ "H have the 
biggest contribution to the source of vector perturbations and hence the magnetic fields. 
Thus we can set g ~ l/r/ in investigating the behavior of the source. When the tight coupling 
approximation is valid, i.e. 2 ; i:|> z^ec, we can estimate the time evolution of 6v as 

5v oc oc r/'^, (4-19) 

T 

from eq. (3.10) and eq. (3.11). Then after recombination, we can estimate 6v ~ from the 
same equations as 

6v oc k^p^a oc (4.20) 


Using these relations, the time evolution of the source term of magnetic fields before and after 
recombination can be derived as 


S{p) = p^Yv^q^^"^ oc 


r/^/^ (before recombination) 
p~^/'^ (after recombination) 


(4.21) 


Therefore, the evolution of magnetic fields becomes source-free after recombination. In fact, 
during recombination, 5v is considerably enhanced and significant amount of magnetic fields 
is produced by the end of recombination 2 ; > 300. For 2 ; < 300, the magnetic fields simply 
decay adiabatically. The spectrum of magnetic fields after recombination is depicted in figure 

3. 
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k[Mpc''] 


Figure 3. Power spectra of magnetic fields at ^ = 2000 (black solid line), 1500 (green solid line), 
1000 (cyan solid line), 700 (green solid line), and 500 (red solid line) from scalar fields in the NLSM. 
Here we take the NLSM parameters as /N = 10“^^mp[. The evolution of the spectrum comes to 
an end by z = 300. 


5 Conclusion 


In this paper, we consider magnetic field generation from self-ordering scalar fields that follow 
the NLSM. We find that to reliably estimate the magnetic fields, one needs to expand the 
Boltzmann equations up to the third order terms in the tight coupling approximation. This is 
because the vorticity a is very large so that O {{k/f)vj) ~ O [{k/f)^a^ in the tight coupling 
era (see figure. 1) if the anisotropic stress of the scalar fields eq. (3.6) is an external source. 
By smoothly connecting the tight coupling solutions to the numerical ones we obtain the 
full magnetic field spectrum in the radiation dominated era and the matter dominated era, 
with an analytic interpretation of the results. In so doing, we see that the scalar fields 
with wavenumber q ^ 1/r] are the main source for both super-horizon {krj <C 1) and sub¬ 
horizon [krj ;§> 1) magnetic fields. By extrapolating our numerical result toward smaller scales 
analytically, we find B ~ 10“^^ (1 -|- (fc/Mpc“^)^^^ /\fN Gauss at A: > 1 

Mpc“^ and z > Zrec- This might serve as a seed of large scale magnetic fields in the present 
universe. 
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A Appendix 

A.l Appendix A : Power spectrum in the NLSM 

Let us consider the power spectrum of some variable X{k) and Y(k) which is generated by 
scalar helds following the NLSM. At first, let us write X{k) using its transfer function Fx{q,p) 
and FY{q,p) as 

^i^) = j (27r^3 i^J^3 ^x{q,p)l3sip)l3s{^{2TTf6{k - p - (A.l) 

and a similar expression for Y{k). The (cross) power spectrum of X and Y is defined by 

(^X{k)Y*{k')) = {2xfVxY{k) 5(k - k'). (A.2) 

To calculate the power spectrum, we need to calculate four point correlation function of scalar 
fields. This is given by 


(/3a(gi)/3a(pi)/3b(g2)/3b(p2)) = (/3a(gi)/3a(pi)) {PbiQ2)l3b{P2)) 

+ Wa{qi)l3b{q2)) {/3a{pi)l3b{p2)) 

+ Wa{qi)l3biP2)) {I3a{pi)/3t{q2)) 

= (Gvr^)^ {2TTfd{pi + qi) (27r)^h(p2 + ^ 2 ) 

(GTT^)^ n q 

+ ^ (2vr)^h(^L - 92 ) {2TTfS{pi -P 2 ) 

(2vr)^h(gi -P 2 ) {2TTf6{pi - 92 )- (A.3) 

Here, = pA is the initial power spectrum of the scalar fields. Substituting eq. (A.3) 

to eqs. (A.l) and (A.2), and using F{q,p) = F{p,q), we calculate the cross correlation as 


X{k)Y*{k')) = 


f d?qi d?pi d?q2 d?P2 

J (27r)^ (27r)3 (27r)3 (27r)^ 

X {27rfd{k - Pi - qi) {2TTf6{k' - P2 - ^2) 
I' (fqi d?pi d^q2 d?‘P2 

J {2tt)^ (27r)^ (27r)^ (27r)^ 

(GTT^)^ q q 

{2TTf6{qi -^ 2 ) (27r)^h(pi - P 2 ) + 


Fxiql,pi)FYiq2,P2) Waiqi)l3aipi)l3biq2)l3biP2)) 

^6{k' -p2- q2) 

Fx{qi,Pi)FY{q2,P2) 


(67r2) 


212 


N ' ^ ^ ^ N 

X {27rf6{k - Pi - qi) {2TTf6{k' - P2 - ^ 2 ) 

= VF{27rf6ik-k') 

f d^q 


{27rf6{qi - P 2 ) {2'nf5{pi - ^ 2 ) 


J (27r)3 
= 2VF{2T,)H{k-k') I 


Fx{q, \k - q])F^{q,\k - q\) + Fx{\k - q],q)FP{q,\k - q]) 
d^q 


(27r)3 

Therefore the power spectrum VxY is read off as 

d^q 


Fx{q,\k - ^)Py(g, \k-^) . 


(A.4) 


VxY{k) = 2V, 


N 

ini 


(2vr) 


Fx{q,\k - ^)Py(g, \k - ^) 




N 


= ^ I dq q'^ dp Fx{q, \k - ^)Py(g ,\k - ^), 


(A.5) 
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where ^ = k ■ q. 
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